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Recent results suggest that the use of ensembles in Statistical Mechanics may not be necessary for 
isolated systems, since typically the states of the Hilbert space would have properties similar to the 
ones of the ensemble. Nevertheless, it is often argued that most of the states of the Hilbert space 
are non-physical and not good descriptions of realistic systems. Therefore, to better understand 
the actual power of typicality it is important to ask if it is also a property of a set of physically 
relevant states. Here we address this issue, studying if and how typicality emerges in the set of 
matrix product states. We show analytically that typicality occurs for the expectation value of 
subsystems' observables when the rank of the matrix product state scales polynomially with the 
size of the system with a power greater than two. We illustrate this result numerically and present 
some indications that typicality may appear already for a linear scaling of the rank of the matrix 
product state. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

Statistical mechanics has been very successful in mak- 
ing predictions about the behavior of macroscopic sys- 
tems we encounter in nature, yet there are still open 
questions concerning a fully quantum formulation of sta- 
tistical mechanics. For instance, how can we explain the 
use of statistical ensembles in the description of a physical 
system which is supposed to be in a definite state? In this 
context the independent re-discovery by several groups 
of the importance of typicality (T-'^ hasgiven rise to in- 
teresting directions of research [5,, ^, |^, |3| ■ Although the 
concept originally appeared as an incomplete formulation 
in a work by Schrodinger (To| . Lebowitz was the one to 
coin the term 'typicality' and, together with others, did 
some pioneering work on the subject 0, Q . Typicality 
can be seen as a key feature justifying the effectiveness 
of standard equilibrium statistical mechanics, without re- 
quiring crgodicity or mixing. The works on typicality in 
the quantum setting have shown that ensemble averages 
and subjective ignorance may not be necessary concepts 
for the understanding of statistical mechanics [3| . 

Intuitively, typicality refers to the fact that the vast 
majority of pure microstates of a quantum system, be- 
longing to a well-defined region of the allowed state space, 
yield measurement outcomes very close to each other. 
More quantitatively, typicality can be associated with a 
very small variance of the measurement outcomes with 
respect to a specified ensemble of states. 

Previous works 0-01 have focused on the study of typi- 
cality for general quantum states, providing a first alter- 
native approach to the foundational problems of quan- 
tum statistical mechanics. However it is well know that 



the generation of Haar distributed random states is hard 
even at the quantum level [71]. Therefore in order to con- 
sider typicality an effective scheme for the justification 
of statistical mechanics one should restrict to realizable 
random states, possibly with some specific physical con- 
tent. First we need to choose and characterize this set of 
states, though of course the choice is not unique. In the 
present work we focus on Matrix Product States (MPS) 
(see [m for a review and original references) as instances 
of physically accessible states. The reason why we re- 
strict the study of typicality to MPS is because these 
are a good example of physically relevant states, arising 
as ground states of local Hamiltonians and being at the 
basis of some of the most recent and promising classical 
algorithms for the simulation of quantum systems [llj. 
Both of these properties justify a better understanding 
of their statistical properties with respect to typicality, 
which eventually can also lead to new powerful simula- 
tion techniques (as the work in Q may suggest). We 
shall prove that typicality can emerge in the MPS set, 
and then illustrate this result with some numerical sim- 
ulations. 



II. RANDOM MATRIX PRODUCT STATES 

A matrix product state is a pure one-dimensional quan- 
tum state whose coefficients are specified by a product of 
matrices. In the case of Periodic Boundary Conditions 
(PBC) an MPS can be written as 



= ^ Tr [A'^ [1] • • • [N]] \ii 
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whereas for Open Boundary Conditions (OBC) we have 
V (0,|A'4l]...A^"[7V]|0f^)|*i...z^), 
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with |(/)/) and \(I)f) specifying the states at the bound- 
aries and life) a local basis at site k. The ma- 
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trices {A^[slA^[s],...,A^[s]}, with s G {l,...,iV}, 
are x-dimensional complex matrices, with D the local 
Hilbert space dimension. For homogeneous MPS the set 
{v4^[s], ^^[s], . . . , A^[s]} is the same for all sites s. In 
the case of PBC they are referred to as Translationally 
Ivariant (TI) MPS. In the present work, for simplicity of 
notation and analysis, we deal numerically with OBC- 
MPS and analytically with PBC-TI-MPS. We checked 
numerically that all of our conclusions hold true indepen- 
dently of the boundary conditions and invariance under 
translations. 

By definition, a MPS is specified by the set 
{A^,A^, . . . , A^}, though there may exist a different set 
of matrices that form the same MPS. In [l^l it is shown 
that this sort of gauge degree of freedom can be fixed us- 
ing a canonical form where the A matrices satisfy two 
constraints: Ez=i ^'^'^ = I and E^Li ^'^^A' = A, 
for fixed A (an alternative set of constraints is given by 

E^i A'^A" = I and A'AA'^ = A, see [H for de- 

tails). MPS can also be seen as generalized valence-bond 
states [llj, and as such, emerging from the projection on 
some virtual or ancillary Hilbert space. The fundamental 
parameter of an MPS is the size x of the A-matrices. In 
general, an MPS contains NDx^ parameters, much less 
than the usual of a general state, and as a conse- 
quence the maximum entanglement a subsystem can have 
with its environment depends on x- It can be shown that 
any state can be described as an MPS for large enough x 
with at most x (though there is no advantage in 

such a description) For more details and properties 
of MPS used in this work see appendix B. 

For our purposes we need to generate an ensemble of 
Random MPS (RMPS) and the way to do this is by no 
means unique. One could think, for example, of choosing 
a set of matrices {A^, . . . , A^} belonging to some relevant 
ensemble known in random matrix theory. This choice 
would induce additional symmetries on the ^-matrices 
that will constrain the set of RMPS too much, and for 
which the physical meaning would not be clear a priori 
(see |12J for a related construction in a different context). 
The ensemble of RMPS that we consider in this work is 
constructed by the repeated random unitary interaction 
between an ancilla and a physical system, as described 
in the framework of the sequential generation of MPS 
[13I [l5| . This is an operationally and physically moti- 
vated realization of MPS. We now briefly summarize the 
construction in [15]. Consider a spin chain initially in a 
product state |0)®^ e "Hf^ (with Hb ^ C") and an an- 
cillary system in the state {(pi) G Ha — 'C^- Let U[k] be 
a unitary operation on Ha ®'Hb, acting on the ancillary 
system and the fc'th site of the chain (see FigH]). The 
A[k] matrices are defined by 



AL,^[/c] = (z,a|C/[/c]|/3,0), 



(1) 



where the greek indeces refer to the ancilla space and 
the latin indeces to the physical space. For homoge- 
neous MPS the unitary interaction is the same for all 
the sites in the spin chain. Due to unitarity we have 



|0> |0> 



|0> 



U[l] 



U[2] 



U[N] 



w> 



Figure 1: Sequential generation of an MPS with \<j)i) and 
\(f)f) the boundary normalized states. 



Y,i A^[k\^A^[k] = for aU k in the bulk. This property, 
together with a proper normalization of the boundaries, 
corresponds to an MPS of unit-norm (see appendix B for 
more details). Letting the ancilla interact sequentially 
with the sites of the chain and assuming that the an- 
cilla decouples in the last step (this can be done without 
loss of generality, as shown in [la|), the state on "H^^ is 
described by 
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which is a homogeneous MPS with OBC. It can be proved 
[13I [l5| that the set of states generated in this way is 
equal to the set of OBC-MPS. We choose the interaction 
characterizing the homogeneous RMPS ensemble to be 
represented by a random unitary matrix U distributed 
according to the Haar measure. 

Since any state can be described by an MPS when 
X oc UMj typicality should appear trivially for MPS 
with this exponential scaling of x in N. Therefore the 
relevant question is if it is possible to have typicality when 
the rank increases at most polynomially with the number 
of particles: x = polyO>i)- This will be the subject of the 
next section. 



III. TYPICALITY IN RMPS 

Typicality can be studied at a more formal level in 
the framework of concentration of measure, a mathemat- 
ical tool which allows to establish typicality in large- 
dimensional Hilbert spaces jlOl]. The concentration of 
measure phenomenon allows to quantify the probability 
of fluctuations for functions of random variables, and in 
the physical literature has already been applied in a vari- 
ety of contexts [l7l - l20j . We shall use a result on the con- 
centration of measure phenomenon for the unitary group 
in order to prove typicality for subsystems' observables 
(a more rigorous mathematical introduction to the topic 
can be found in (l6j : in particular, here we use theorem 
6.7.1 of that book). Concentration of measure holds for 
the unitary group and this means that there exist univer- 
sal positive constants ci and C2 such that for any function 
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/ : U{d) — >■ M, from the set of Haar-distributed unitary 
matrices of size dx d into R, and with Lipschitz constant 
77 



Pr [|/ - f\ > e] < ci exp {-C2e^d/r]'') 



where / denotes the average value of /. From this ex- 
pression one can see that typicahty is vahd only for func- 
tions / for which the ratio d/rj'^ (where 77 can in principle 
also depend on d) increases with the dimension d of the 
domain, since in this case the probability of large fluc- 
tuations around their average will decrease exponentially 
in d. 

In the present work the random variable / will be the 
expectation value of an observable with respect to a ran- 
dom MPS. The observables that we consider are those 
that can be expressed as the tensor product of local ob- 
servables. In the usual transfer matrix notation for nor- 
malized MPS (see appendix B) we can write 



/ = Tr 



N 



.k=l 



0[k] 



with 



D 



Eo[k]= J2 {^k\0[k]\Jk)A'nk]®A^nkr 

ik,jk = l 

the tranfer matrix associated to the observable O = 
^k=i ^i^]- 74[fc]-matrices characterizing the state 

are obtained as sub-blocks of random unitaries U[k], ana- 
logusly to Eq.([T]). In this way the expectation value of 
the tensor product of local observables can be seen as a 
random variable / : U{xD) — >■ R, from the set of uni- 
formly distributed unitary matrices of size xD x xD into 
M. In order to apply the concentration of measure re- 
sult for functions of random unitaries we need to find an 
upper-bound for the Lipschitz constant 77 in ([2]) 



77 = supu^^u^ 



\f {Ui)-f{U2)\ 
\U1-U2h 



In order to do that we look for an upper-bound of the 
absolute value of the differential of / 



\df\ = \dTr 



N 

n ^oik] 



.k=l 



We consider the case of subsystems of size L specified by 
observables O of this form 



N 



o = 



)0[k] 



I[k] 



differential calculus for matrices we have (see appendix 
A) 



(2) \df\ 



dTr[[[Eo[k]E, 



N-Li 



k=l 



< 



Tr 



df[Eo[k])E, 



k=l 



Tr 



W Eo[k]dEj 



N-L 



.k=l 



(3) 



Let us consider the first term on the right-hand side in 
We have that 



Tr 



dY[E, 

\ k=i 



Oik] I Ej 



N-L 



L L 

<J2I[\\^o[j]\\oo\\dEo[k]\U\E, 

k=l j^k 



N-L] 



where \\-\\k stands for the matrix fc-norm. We denote with 
||0||oo the max {||0[fc] ||oo ; k = I, . . . , L} and use the fol- 
lowing relations 

Ploo = ||(K)(i|s®lA)C/(|0)(0|B ®Ia)||oo < \\U\\oo = 1, 



\E 



N-Li 



< 1 + e 



2 N-L 
2 I 



\\dEo[j]\\oo < 2D^\\0\\a,\\dU\\oo, 

where £2 is the second-largest eigenvalue in the spectrum 
of the transfer matrix over the ensemble of RMPS (re- 
member that for normalized MPS the largest eigenvalue 
of -Ei is 1, so that €2 < I for all the realizations). We can 
then bound the first term of ^ with 

2LD'^{l + x'e^~'^)\\0\a\dU\U (4) 

Let us now consider the second term in the right hand 
side of ([3]), which can be written as 



N 

E 

3=L+1 



Tr 



.k=l 



Kk=L+l 



We can bound each term in this sum with 

L _ 

II II i?owi?|^^di?i[,]||oo|l4^^||l, 

and the total sum with 



We will refer to the rest of the chain, of N—L sites, as the 
environment or bath. Using standard properties of the 



/ + X'4 ^ ^){N~L~l)\\dE,\\^\\Y[Eo[k]\\oo, 
^ fc=i 
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which is smaller than or equal to 

2{N ~L - 1)D^^+'^ 



Empirical probability distribution 



I N-L I 

1 + x' 



lOll^lldL/ll, 



(5) 

The total variation of the expectation value is then 
bounded by the sum of equation (jH) and equation ([5]), 
which in the relevant regime of interest {N ^ L) be- 
comes 



2L+2 



^l|0||^||dC/||2, 



where we used the fact that ||-||oo<||'||2- The Lipschitz 
constant rj for the function / is then upper-bounded by 



Tj < 4D^^+^N\\0\ 
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for ^ L. Along with Eq.®, this implies that increas- 
ing the size of the environment will cause the expectation 
value of the observables of any subsystem to concentrate, 
provided x(A^) increases faster than N'^ 

Pr [|/ - f\ > e] < ci exp {~c!,e\{N)/N') , 

where we have absorbed all the constant in It is im- 
portant to notice that the set of MPS, for fixed x ^iid 
N, is exponentially small with respect to the total num- 
ber of states in the same Hilbert space. As N increases 
the dimension of the Hilbert space will increase exponen- 
tially but a polynomial scaling of x in -/V will be sufficient 
to guarantee typicality. This shows that typicality is a 
property of a class of accessible states of quantum system, 
extending some of the implications of previous work [3- 
Q to an experimentally and computationally accessible 
regime. 

Trace distance from mixed state 




Figure 2: Trace distance between the average OBC-MPS with 
X = 2 and the completely mixed state. The lower line is 
the trace distance between the average uniformely distributed 
random states and the completely mixed state. 



IV. NUMERICAL RESULTS 



The recent new approach to typicalit; 
main flavors, one due to Popescu et. al. 



comes in two 
and a second 



Figure 3: Empirical probability distribution for the eigenval- 
ues of the reduced density matrix of RMPS (5000 states). The 
subsystem has dimension 2 and the total system has dimen- 
sion 16 (x = 2). 



one due to Reimann [3| . The first approach is more math- 
ematical in two aspects: it uses results from the concen- 
tration of measure phenomenon and considers distances 
between states of subsystems. Reimann, on the other 
hand, uses more heuristic arguments and studies general 
observables. In this latter approach one has to show that 
the variance of the expectation values of observables is 
small and decreasing with the size of the bath. This 
implies, by the Chebyshev inequality, a concentration of 
measure result. With the approach of [4] one can directly 
study the fluctuations of the trace distance between the 
states in the ensemble and their average. A concentra- 
tion result with this approach is a stronger result, in the 
sense that it is sufficient for having typicality at the level 
of observables. On the other hand typicality for all local 
observables can also imply a weaker concentration result 
for the state of the subsystem (see appendix C). In the 
numerical simulations we considered both the variance of 
the expectation value of observables and the fluctuations 
in the subsyste m of th e distance of RMPS from their av- 
erage state 75 = for which an exact expression can 
be obtained and reads in terms of the components 



Pi i — Tr.ij®2iv 



(|0, j) (i, 0\(E)T(E) T^) C/®^ (g) C/t®^ 



where the A^— components vectors are defined by = 
(0, 0), i = (ii, . . . , iat) and the same for j. The operator 
T acts on the A^— tensor product of the ancilla system and 
cyclicly permutes the components 



T\a) 



\a)N = \a)N\a)i . . . \a)N- 



The unitary random matrix U{xD) is the one used in 
the sequential generation of the random MPS. A closed 
form for the average of the tensor product of unitaries 
U»N ^ [/t«.JV jg known [IJ. In Fig (H we plot the trace 
distance between the average random MPS with OBC 
and the completely mixed state (for a chain of 4 qubits) 
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as a function of the size of the samphng set 
(the number of randomly generated states). In the same 
figure we also plot the same quantity in the case of ran- 
dom general pure state (not necessarily MPS). As can 
be seen from the figure the average OBC-RMPS is at 
a finite distance from the mixed state, while the aver- 
age general state approaches the mixed state increasing 
the number of sampled states. Another distinctive fea- 
ture of the homogeneous OBC-RMPS is shown in Fig|31 
where we plot the empirical probability distribution of 
the eigenvalues of the reduced density matrix of 1 qubit 
in a 4-sites RMPS. In Fig|3]we show the same plot for 
general randomly generated states. As can be seen the 
two distribution differs significantly. From now on, un- 
less otherwise stated, for all simulations we consider an 
ensemble of 500 RMPS which originate from random uni- 
taries distributed according to the Haar measure. We 
now want to illustrate our analytical results studying the 
behavior of the variance of a^, which acts on a particle 
(i = 1) in the middle of the chain. Note that the im- 
portant variable is not the absolute size of the subsystem 
or the bath but the ratio between them. As can be seen 
in Figinia, when we fix the value of x and increase the 
number of qubits in the bath the variance starts to de- 
crease, but soon reaches a limiting value and does not 
decrease any more. The limiting value depends on 
becoming smaller as x is increased. This could be ex- 
pected from our bound and from a known property of 
MPS: correlations between system and environment are 
of finite range and depend on x- This is also consisten 
with recent result on finite entanglement scaling at criti- 
cality [22|. Note, however, that our analytical result does 
not exclude the possibility of having typicality for fixed 
X 01 x scaling linearly with N. It only guarantees typ- 
icality in the case of a scaling greater than quadratic. 
We then analyze the case where x — ^ L, as show in 
FiglSlb. There it can be seen that until x — N — L — 180 
the variance is decreasing monotonically, which indicates 
that typicality can emerge already for a linear scaling of 
X with the number of particles. However, at the present 
moment our simulations do not allow for a conclusive 
statement about the precise scaling of x with N — L that 
assures typicality. 

We also investigate the behavior of the average trace 
distance from the average state at the level of the sub- 
system, which is one particle in our case. Denoting with 
Ps the reduced density matrix of the subsystem of a 
RMPS and pj the average MPS, Fig. Pa) shows the de- 
pendence on the size of the bath of the average value 
E{Di) of the trace distance Di = \\ps — plUi. In general, 
we expect that having typicality at the level of states 
is harder than at the level of observables (see appendix 
C). Again we look at the case of fixed x (Fig-(ISla)) and 
X — N — L (Fig.(inib)). The conclusions are similar to 
the case of observables, and it appears that even at the 
level of states typicality may occur already for a linear 
scaling of x in the size of the bath, while the behavior 
for different but fixed values of x is again consistent with 



Empirical probability distributipn 
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Figure 4: Empirical probability distribution for the eigen- 
values of the reduced density matrix of a general Haar- 
distributed random state (5000 reahzations). The subsystem 
has dimension 2 and the total system has dimension 16. 




Figure 5: (a) The variance of the expectation value of {L = 
1) increasing the size of the system and for fixed but different 
values of X = 2, 4, 6, 8, 12, 20, 50, 100, and 180 (from top to 
bottom), (b) The variance of the expectation value of (Jx 
[L = 1) for increasing system size when the MPS dimension 
increases linearly with the number of particles in the bath: 
X = N-L. 



well-known MPS properties. 



V. CONCLUSIONS 

In summary, we have shown that typicality can arise 
not only for an exponentially big Hilbert space but also 
for a physically accessible smaller set of states: the Ma- 
trix Product States. More specifically, we showed ana- 
lytically that typicality occurs for MPSs having x scal- 
ing polynomially in the size of the system (with a power 
greater than 2). We then presented some numerical 
calculations which indicate that typicality may already 
emerge for a linear scaling of x in the system size at the 
level of both observables and the states. 

Our results provide further evidence that typicality 
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Figure 6: (a) The average trace distance from the average 
state of the ensemble of RMPS as the number of qubits in 
the environment (A'' — L) increases, for fixed but different 
values of X = 2, 4, 6, 8, 12, 20, 50, 100, and 200 (from top to 
bottom) and for L = 1. (b) The average trace distance from 
the average state of the ensemble of RMPS as the number of 
qubits in the environment (A'' — L) increases, taking the MPS 
dimension to increase linearly with the number of particles in 
the bath: x = N ~ L. 



may play a role in a better understanding of the foun- 
dations of statistical mechanics. Nonetheless, there are 
still some aspects that require a deeper analysis. For ex- 
ample, it would be interesting to have more information 
on the average state obtained from the present ensemble 
of RMPS (in general p will still be a matrix product state 
but with a bigger rank than its components) and to bet- 
ter characterize the role of the geometry of the partition 
in between subsystem and bath and their correlations. 
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APPENDIX A. If / is a real-valued function oir a 
metric space {X, d) , its Lipschitz constant is 



\f\c = SUp^^y 



d{x,y) 

for x,y € X. If X = R", we let d{x, y) ^ \\x - y\\, with 




the euclidean norm. The above notation can be easily 
applied also to the case when X — A^„(C) is the set of 
n X n matrices. For A e A^„(C) we define 



I^IV^)II 



lAh^^TriAfA), 



\\Ah=Tr\A\. 

All this norms are unitarily invariant and submultiplica- 
tive. In this work we used also the following relation 

Tr{AB) < ||A||eo||S||i. 
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In the derivation of the upper bound for the Lipschitz 
constant we made use of standard properties of differen- 
tiation with respect to a matrix (|23|) 

dTr{X) = Tr{dX), 



d{XY) = {dX)Y X{dY), 

where X and Y are arbitrary matrices. Defining the dif- 
ferential df{X) to be the part of f{X + dX)-f{X) which 
is hnear in dX, the gradient V/ satisfies df = V/ • dX. 
Since the Lipschitz constant is equivalent to the sup|V/|, 
an upper-bound for the gradient will provide an upper- 
bound for the Lipschitz constant. 

APPENDIX B. Here we review some of the notation 
used in the literature on MPS. A detailed exposition can 
be found in [ll| . For brevity we will focus for the moment 
on normalized MPS with periodic boundary conditions, 
which by definition is a state that can be written as 

D 

IV) = Tr{A^'m'm---A^"[N])\iit2...tN), 

ii,...,ijv = l 

where the matrices A are x^X complex matrices, labeled 
by the site index e {1, ■ ■ ■ ,N} and by the local bases 
index € {1, . . .,£>}. 

The expectation value of some operator S, which is the 
tensor product of local operators S[k] at each site k 

S ^ S[1]<S)---®S[N], 
is given by (■01 (^^j^ S'[/c]|'0) which is equal to 

f2 Tr(Y[A^''[k]\Tr(f[A^Hkr) 

ii,i^...,2iv,i^ \k—l / \k—l / 



N 



k=l 



N 



N 



J2 TriJlA'" [k] (g, Y[ J^'" [k]* 

\k—l 



Tr 



N D 



\{Y.^'k\Smk) [A^nk]®AHk\ 

fc=l ik,i't. 

Defining the transfer matrix or transfer operator 

D 

Es^^m ^ Y.^'k\S[k]\h) [A^''[k]®A^'>'[kX 
we see that 



N 



)S[k]\^)=Tr 



fe=i 



AT 



l[^s[k][k] 



The normalization of the state is given by 

JV 



Trill E^,][k]), 



(6) 



fe=i 



fe=i 



with S[k] = l[k] for all k. 

Let us now consider the case of sequentially generated 
OBC-MPS [1^ . Consider a spin chain initially in a prod- 
uct state 

and an ancillary system initially in the state "Ha- 
Let us introduce a unitary operator U[k] acting on Ha ^ 
Hb, for each site k in the chain. Defining 

Al,^p[k]^{i,a\U[k]\(3,0), 

unitarity implies the following 

C/[fc]^C/[fc] = ^A*[fc]U'[/c] =1;, ke{l,...,N}. 

i 

(7) 

Letting the ancilla interact sequentially with all the sites 
in the chain (see FigH]) we obtain the state 

D 

ii,...,iiv=l 

Let us write the OBC-MPS in a different way (A*"+i [N+ 
l] = ((/.f|and A'o[O] = |0,)) 



JV 



\{'^'k\smk) 



k=l 



D 



ll,lj...,ljV,ljv 



N 

Y[ [fc] (g)yl*4fc]' 



Tr 



N 



l[('''^\S[k]\tk) 



k=l 



D N 

E l[(''k\S[k]\tk){A^nk]®A^'''[kr 



i,i{...,iN ,i'f^ k=l 



jl,...,j]V = l 



xA'^iNjA'^miN+iiN-.-iiio), 

where A''' [k] are D^+i x Dk matrices with Dq = D^+i = 
1 and Dk = X with k G {1, . . . , N}. The normalization 
of the MPS is given by 



(VIV) 



D 

E 



A'" [0]U*4l]^..^'"+i[iV+l]1'x 
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X [N + l]A'" [N]... A'^ [1]A'« [0]. 

Using a singular value decomposition [isj it is always 
possible to find a canonical form for !■(/;) such that 

[N + [iV + 1] = I-^. 

Using ([7]) recursively one has 

(VIV') = Y,A''>[0]''A'°[0] = 1, 

io 

where the last equality follows from a normalization con- 
dition which can be imposed without loss of generality 
on the boundary local matrix A* [0] . 

In the case of periodic boundary conditions the norm 
of the state 



J2 TriA''[l]A''[2]...A"'[N])\i 



ii2 ■ ■ -In), 



ii,...,ijv = l 



is given by EqlH with Ei[k] = E^=i^"=W A^^ik]*. 
For simplicity of notation, but without loss of generality, 
we shall restrict now to the translational invariant case. 
To any MPS it can always be associated a Completely 
Positive (CP) map [H 



£{x) = J2a^xa. 



The CP map can always be assumed to have spectral 
radius 1 (corresponding to the absolute value of its max- 
imum singular value) [13|]. The map £ and Ei have the 
same spectrum since 

(/3i|£(|ai)(a2|)|/32) = /^zl^ilai, aa). 

Assuming that £ has only one eigenvalue equal to 1 (with- 
out loss of generality [l^ ) , one see that for N big enough 



where Aiis the maximum eigenvalue equal to 1. The cor- 
rections are exponentially suppressed in N. In our an- 
alytic derivation of an upper bound for the Lypschtiz 
constant of the function 



we assume TV to be big enough that for all purposes 

N 



Trill Esik][k]), 



k=l 



without the need of any additional normalization (and 
numerically this is verified to be true already for iV > 4) . 



APPENDIX C. A concentration of measure result 
obtained for the trace distance between random states 
and their average implies the same concentration of mea- 
sure result for the expectation values. This can be proved 
using the following general relation between the trace dis- 
tance of normalized states and the difference between the 
expectation values of an observable A 

\TripA)-Tr{M)\ = \Trip~^)A\ < ||p - VlliPlloo, 

and we assume the operator norm of A finite. If p is a 
random state and i/; is its average, one can see that a 
bound on the fluctuations of the right hand side of the 
above inequality implies a bound on the fluctuations of 
the left hand side. In general one can say that close states 
will have close expectation values for any observable of 
finite operator norm. But in general the converse is not 
true: if the expectation value of some observable with 
respect to different states is close this does not imply that 
the states are close. A way to estimate the state-distance 
is shown in [4]. 

We prove in this manuscript that for the expectation 
value / of any local observable, restricted to a subsystem 
of size L much smaller than the size N of the total system, 
the following holds true 

^^[1/ -71 > e] < ciexpi-c'^e^/N^), 

where the sampling is done with respect to a set of ran- 
dom matrix product states of rank x- Without lack of 
generality lets restrict to a chain of qubits. Any operator 
in the subsystem of size L can be expressed in a basis 
of 4^ unitary orthogonal operators (for an explicit con- 
struction see [1|). Lets call each element in this basis Ux, 
with xe{l,...,4^}. We shall indicate with p a realiza- 
tion of a normalized random matrix product state and 
with p the average state. Lets define Px{p) = Tr{Uxp) 
and Px(j}) = Tr[Ux'p)- The previous concentration result 
holds true in particular for px (for any x) 

Pr[\px{p) -Px{'p)\ > e] < ciexpi-d^e^x/N^)- 

Since there are 4^ different x values we can also writhe 

Pr[3x : \Px{p) - Px{-p)\ > e] < 4^ciexp(-4e\/iV2). 

Any random density matrix associated to the normalized 
RMPS can be written as p — ^xPx{p)Ux- When \px{p) — 
Pa;(p)| < e for all x, then it follows 



IP- P\\2 



YiPxip) - Px{p))U, 



x\\2 



H 2 



X! ^Px{.P) - Px{p)) Ux 



=4^5^(p.(p)-p.(p))2<42V 

X 

We can write then 

||p-p||i<V4^||p-p||2< 43^/2^ 
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and from this it follows 

Pr[\\p-p\\i > 4'''/'e] < 4^ciea;p(-4e\/Ar2). 

Which tells us that keeping L fixed and much smaller 
than N, for x oc AT^and e oc A''~^/^one has 



This prove a concentration of measure result "at the level 
of states". But this result is weaker with respect to the 
result for the observables by a factor 4^. 



